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£NJ ■ We study the eigenvalue problem of the rational Calogero model with the coupling of the inverse- 

square interaction as a complex number. We show that although this model is manifestly non- 



invariant under the combined parity and time-reversal symmetry VT, the eigenstates corresponding 
to the zero value of the generalized angular momentum have real energies. 
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The standard practice in quantum mechanics is to consider self-adjoint operators so that 
the corresponding spectrum is real and consequently, the time-evolution of the states is 
unitary. However, it has been found recently that a non-Hermitian Hamiltonian invariant 
under the combined VT symmetry can have real spectrum . The 

spectrum is entirely real if the VT symmetry is unbroken. On the other hand, part of the 
spectrum is real if the VT symmetry is broken spontaneously [l ] . 

In this Letter we study the rational Calogero model Q,Q2| with complex coupling for 
the many-body interaction. Several extensions of the Calogero model with non-Hermitian 
Hamiltonians have been studied in the literature 0, Q] . However, all these models respect 
VT symmetry and lead to real spectrum. The Hamiltonian we study is identical to that of 
the original iV-body rational Calogero model except that the coefficient of the inverse square 
interaction is complex. This model is non-Hermitian, manifestly non-invariant under VT , 
but nevertheless is shown to admit a partly real spectrum, with the ground state energy 
always real. 

Following the analysis of Calogero the eigenvalue equation of the rational Calogero 
model with complex coupling for the many-body inverse-square interaction can be reduced 
to the eigenvalue equation of an effective single particle Hamiltonian containing harmonic 
and inverse- square interaction. We show that for certain values of the complex coupling 
of the inverse square interaction, this effective Hamiltonian can indeed be Hermitian when 
the generalized angular momentum ^j] is zero. This however happens only for the strongly 
attractive values of the inverse square coupling term in the effective Hamiltonian. Such a 
strongly attractive single particle system in absence of the harmonic term has already been 
analyzed by Case Q], who obtained the corresponding spectrum in terms of an undeter- 
mined parameter which has the physical interpretation of a cutoff in the coordinate space. 
In our analysis, we explicitly introduce a cutoff which prevents the particles to occupy the 
same position simultaneously, i.e. it serves as a cutoff in the coordinate space. In the anal- 

n 

ysis of ref. [13], the spectrum of the strongly attractive model in absence of the harmonic 
term is unbounded from below. In the model under consideration here, which includes the 
harmonic term, the spectrum is however bounded from below for finite nonzero values of 
the cutoff. 
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The Hamiltonian of the rational Calogero model is given by 
H = 



N Q2 



i=l ^ 



;X\; 



(1) 



where a and f2 are constants, Xi is the coordinate of the z th particle and units have been 
chosen such that 2mh~ 2 = 1. The constant a is conventionally taken to be real, which is 
required for the operator H to be Hermitian. In this Letter we shall however show that for 
certain complex values of a, the eigenvalue problem 



(2) 



still admits a real spectrum with normalizable solutions. In the discussion below, we there- 
fore take a = or-Mot where or and ai are the real and imaginary parts of a. The Hamiltonian 
H is invariant under parity Xi — > —Xi, while non-invariant under the time- reversal symmetry 
i — > —i. Thus, the combined VT symmetry is not respected by H. 



Following we consider the above eigenvalue equation in a sector of configuration 
space corresponding to a definite ordering of particles given by X\ > X2 > ■ • • > x^. The 
translation-invariant eigenfunctions of the Hamiltonian H can be written as 



* = n 0** - 



(r) P k (x) 



i<j 



where x = (xi,x%, . . . , xn), 



1 



AT ^2( Xi X j) 



N 



(3) 



(4) 



i<j 



and Pk(x) is a translation-invariant as well as homogeneous polynomial of degree k(> 0) 
which satisfies the equation 



f d 2 | ^ 2(q+l) d 

i=l iytj \ X l ~ X j) dXi 



Pk(x) = 0. 



(5) 



The existence of complete solutions of © for real a has been discussed by Calogero . For 
both real and complex a, Pq(x) is a constant and is a solution of Eqn. (5). For k ^ 0, the 
solutions of Eqn. JSJ) for complex a are obtained by analytic continuation in the parameter 
space a. 

Substituting Eqn. J3J) in Eqn. (J2J), using Eqns. (j3J-|2J) and making a further substitution 

_ r -(^+2)$ we get 

H<& = £$, (6) 



where 

Tj u . Kh 1 a) 2 2 /wN 



with u; 2 = ±fi 2 iV and 



H = k + ~(iV - 3) + ~iV(iV - 1)K + ia 7 + i). (8) 

The operator if can be interpreted as the effective Hamiltonian of a particle in the combined 
harmonic plus inverse-square interaction. Let us now choose the real part of a as 

iV-3 1 

flH " ~N{N-1) ~ 2' (9) 



and keep a/ as an arbitrary r ea 
integrable for an + \ > - 



number. It may be noted that the wavefunction \I/ is square- 



ly |1J| , which is indeed satisfied for the above choice of clr for 
N > 3. With the above choice of the complex coupling a, for k = we see that 

fM=^N(N-l)ai = iv (10) 

is purely imaginary and the operator H takes the form 

d 2 iy 1 + 1) 

which is Hermitian and is expected to have real eigenvalues. However, for the same choice 
of or, when k ^ 0, jj, is in general a complex quantity. Thus the eigenvalues of H in the 
k ^ sector in general would not be real. Below we shall first find the eigenvalues of the 
operator H for k = and then provide a suitable interpretation of the states in the k ^ 
sector. 

We now proceed to solve the eigenvalue problem given in Eqn. © with H given by Eqn. 
(jlljl . In this case, the inverse- square interaction is necessarily in the strongly attractive 
regime for v ^ 0. As mentioned before, we put a cutoff r in the coordinate space which 
leads to the boundary condition $(r) = for r = r . We also demand that $(r) e L 2 (i? + ). 
Defining g = o;r 2 and $ = q~^x{q)-> the eigenvalue equation HQ = can be written as, 



d 2 X 



+ 



- + — (l + v 2 ) + — 

4 4g 2 4u;g 



X = 0, X (qo = wr 2 Q ) = 0. (12) 



dq 2 

The above equation can be identified as the Whittaker's equation, the two linearly indepen- 
dent solutions of which are given by 

WE_, ±f (q) = e-*q— M —,l±iu,q), (13) 



where M f^4p — 1 ± iv, q) is Kummer's function jlfij . The general solution of Eqn. (|12|) 
satisfying the boundary condition at go is given by 

X(q) = A \\Ve_ Uq) We_ __-(g ) - We_ We_ ^(g )l , (14) 



where A is the normalization constant. In the limit q — » oo, we have 
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M(a, 6, g) _ e-V a + eV^ + Q (\\ ^ 



T(b) T(b-a) r(a) \g / 



Thus, as g — ► oo, 

X(g) -> Ae2 



r ( 1 + zz/ ) w / n r(i-«/) 

_^(g ) - wj ^_i, ^(go) 



r, 



(16) 



The square integrability of \ can thus be ensured if the quantity in the parenthesis in Eqn. 
(16) is identically zero, i.e. 

Wz%(qa) r(i-«/)r(^- ! 



WjLizfqo) T(l + iu)T( 

A . . 1 '1 V ' V ' * 



-2 1^ / 17 \ 



Aui > 2 



4w ' 



We are interested in the solution of the eigenvalue equation in t 



re limit of a small cutoff, 



l-iv f 

i.e. when g — > 0. In this limit, we have We_ _^(go) — > g 2 (l6jl. The energy eigenvalues 
are thus determined from the equation, 

- iv _ r(i - fr) r(^ - g) 

"° -r(i + w)r(^-|;)' (18) 

or equivalently from 

where and a are the arguments in the polar representation of T(l — iv) and r(iy^ — j-) 
respectively. The spectrum is obtained by solving Eqn. (19) graphically. We have studied 
Eqn. (19) using Mathematica for different values of v and go- For a fixed v and go, there is 
one negative energy bound state and infinitely many positive energy bound states. Moreover, 
the energy levels are not equispaced. The absolute value of the negative energy depends on 
the choice of go- For smaller values for g , \E\ for the negative energy eigenstate increases. 
We plot our results in Fig.l and Fig. 2. 

The following points about the spectrum in the k = sector may be noted: 
1) As shown above, even for certain complex values of the inverse-square coupling, the 
rational Calogero model admits a real spectrum in the k = sector. The real part of the 
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FIG. 1: A plot of Eqn. (JUJ) with w = 0.25, 
v = 0.3 and go = 0.1. The horizontal 
straight line corresponds to the value of the 
l.h.s of Eqn. fgft and / = 2a. 



FIG. 2: A plot of Eqn. (Q]5J) with w = 0.25, 
z/ = 0.3 and go = .03. The negative energy 
state has a lower eigenvalue due to the choice 
of a smaller value of the cutoff go- 



relevant coupling is N dependent and is given by Eqn. (9) whereas the imaginary part of 
the coupling is arbitrary. Thus there is a range of values of the parameter a for which the 
spectrum in the k = sector is real. 

2) The spectrum here for a given value of go and v is bounded from below. In order to see 
this, consider E = —E with S — > oo. In this limit, the argument of the T(^ E - — Jj) can be 
approximated by + ■§-) and Eqn. (18) can be written as 

- v lng - 26 = Q + |-) , (20) 

where ip denotes the digamma function ilfjj. For positive values of its argument, the digamma 
function is monotonically increasing. Thus, for any given finite but nonzero value of go and 
v, E — >• — oo is not a part of the spectrum. 

3) This system generically admits a single negative energy eigenstate. Existence of negative 
energy eigenstates in the Calogero model is associated with the self-adjoint extension of the 
corresponding Hamiltonian Q,^!. We have however not studied the self-adjointness of the 
effective Hamiltonian H here. 

As mentioned before, the parameter /i is complex when k ^ 0. In that case the effective 
Hamiltonian would not be Hermitian and the states with k ^ are expected to have 
complex eigenvalues. These states with complex values of the energy would decay under time 
evolution. Thus the spectrum obtained here has a part which is stable under time evolution 



and a decaying part corresponding to complex eigenvalues. This situation is similar to what 
happens in the case of many nuclei, where only a few energy states are stable and the others 
are resonances with characteristic decay widths. It would be interesting to find such physical 
applications of the model described here. 
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